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Abstract. We study the characteristic size and shape of idealized blazar-induced cascade halos 
in the 1 — 100 GeV energy range assuming various non-helical and helical configurations for the 
intergalactic magnetic field (IGME). While the magnetic field creates an extended halo, the helicity 
provides the halo with a twist. Under simplifying assumptions, we assess the parameter regimes 
for which it is possible to measure the size and shape of the halo from a single source and then to 
deduce properties of the IGME. We find that blazar halo measurements with an experiment similar to 
Eermi-LAT are best suited to probe a helical magnetic field with strength and coherence length today 
in the ranges 10“^^ < So/Gauss < 10“^^ and lOMpc ^ A < lOGpc where % ^ Sq/A is the 
magnetic helicity density. Stronger magnetic fields or smaller coherence scales can still potentially 
be investigated, but the connection between the halo morphology and the magnetic field properties is 
more involved. Weaker magnetic fields or longer coherence scales require high photon statistics or 
superior angular resolution. 
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1 Introduction 

Several lines of reasoning suggest that the voids between galaxies and galaxy clusters contain a 
large-scale intergalactic magnetic field (IGMF). The presence of even a weak IGMF is sufficient 
to explain the observed micro-Gauss galactic and cluster magnetic fields [1-3] through the dynamo 
amplification. Moreover, the IGMF may be a remnant from the early universe as it could have been 
generated during cosmological phase transitions [4, 5], in the epoch of matter-genesis [6-8], or in 
certain inflationary scenarios [9]. An observation of the IGMF and measurement of its energy and 
helicity spectra, therefore, could serve as a powerful new probe of astrophysics, particle physics, and 
early universe cosmology. (For recent reviews on cosmic magnetic fields, see Ref. [10, 11].) 

On the observational side, TeV blazars offer one of the best strategies for measurements of 
the IGMF at redshifts out io z ^ 1 [12-16]. The blazar initiates an electromagnetic cascade as 
its TeV gamma rays produce electron and positron pairs (leptons) upon scattering on extragalactic 
background light (EBL). The cascade develops as the leptons inverse-Compton scatter on cosmic 
microwave background (CMB) photons producing secondary GeV gamma rays. Since the initial 
TeV gamma ray has a mean free path of > 10 Mpc, the cascade develops outside of the host halo 
where it probes the IGMF. In the presence of an IGMF, the charged leptons are deflected by the 
magnetic field and the blazar acquires a halo of GeV gamma rays. In the weak field regime, the 
leptons experience a gentle deflection and the angular extent of the magnetically broadened cascade 
(MBC) goes as © oc S. In the strong field regime, the leptons are so dramatically deflected that the 
secondary emission is isotropized into an extended pair-halo (PH). 

The hunt for cascade halos is ongoing at a number of gamma ray observatories [17-19] as 
well as independent collaborations [20-24]. Most recently Ref. [24] reports evidence for GeV halos 
around 24 low redshift blazars, which are revealed in a stacked analysis of the Fermi-LAT gamma 
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ray data. Assuming a magnetically broadened cascade, i.e. weak bending approximation, they infer 
the IGMF field strength to be So ^ 10“^^ — 10“^^ G where So is the magnetic field strength at the 
present cosmological epoch. 

While there has been extensive analytic and numerical work on the relationship between halo 
size and magnetic field strength, the literature contains little discussion of what information could 
be extracted from the halo shape. In principle the halo size and shape together (morphology) may 
encode not only the magnetic field strength but also its helicity^ 

There are two compelling reasons for considering a helical IGMF. First, magnetic helicity is a 
prediction of many models of magnetogenesis from the matter-genesis epoch [6-8], cosmological in¬ 
flation [25-30], and other early universe scenarios [31-33]. Second, the evolution of a magnetic field 
within the magneto-hydrodynamic (MHD) approximation is known to be more robust to dissipation 
if it is helical [34]. So, if a magnetic field was generated by causal processes in the early universe, it 
has a much better chance of surviving if it is helical. 

The helicity of the IGMF can be observed if charged particles propagate in this field, such 
as in the case of cosmic rays [35] and cascade gamma rays [36, 37]. The parity violating IGMF 
helicity leads to certain non-trivial parity-odd correlation functions of the arrival directions of cosmic 
rays and cascade gamma rays. An advantage of seeking a parity odd signature is that other sources 
of noise are expected to be parity even and hence do not contribute to the signature on average. The 
parity-odd correlation, called Q, has been evaluated for diffuse gamma ray data obtained by the Fermi 
Gamma Ray Telescope, and provides evidence for an IGMF of strength ^ 10“^^ G on distance scales 
^ 10 Mpc with left-handed helicity [38, 39]. 

In this work we endeavor to develop an understanding of how a helical IGMF leads to parity- 
violating features in the shape of blazar-induced cascade halos. To that end, we study an idealized 
system: we focus on a few toy models of the IGMF, we do not model astrophysical sources, we do not 
include stochasticity in the development of the cascade, and we do not include foreground (noise). 
In this setting we can focus on the relationships between the parameters of the magnetic field model, 
Le. the field strength, coherence length, and helicity, and the resulting size and shape of the cascade 
halo. We find that there are several regimes in the range of interesting physical parameters that each 
lead to qualitatively different halo morphology. We identify the region of parameter space where 
measurements of cascade halo size and shape are best suited to probe the helical IGMF. 

In Sec. 2 we review the physics giving rise to the cascade halo. In Sec. 3 we study the prop¬ 
agation of particles through the cascade and derive a set of equations that can be solved to find the 
size and shape (morphology) of the halo. In Sec. 4 we apply the results of Sec. 3 to calculate the halo 
shape for five specific magnetic field configurations, three non-helical and two helical. In Sec. 5 we 
analyze which regions of the magnetic field parameter space could be probed by measurements of 
cascade halos. We conclude in Sec. 6 with a discussion of potential directions for future work. 

2 The Blazar-induced Cascade Halo 

As depicted in Fig. 1, a blazar’s TeV gamma rays initiate an electromagnetic cascade when they 
scatter on extragalactic background light. If the cascade occurs in the presence of a magnetic field, 
the blazar acquires a halo of GeV photons. Then the shape and angular extent of the halo are related 
to the magnetic field in the neighborhood of the blazar. This section is a review of the physics giving 
rise to the halo following Ref. [40]. 

^ A helical magnetic field has a larger amplitude in either left- or right-circularly polarized modes. 
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Figure 1. The geometry of gamma ray propagation. A TeV gamma ray (red arrow) travels off of the line of 
sight (black dashed line) connecting the blazar (yellow star) to Earth (pale blue dot). Pair production occurs 
and either the electron or positron (green arrow) is deflected back toward the line of sight. Inverse Compton 
scattering creates a GeV photon (blue arrow) that eventually reaches Earth. The lepton need not make a 
complete orbit as we have shown here. 


The comoving distance from the Earth to a blazar (source) at redshift is 
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where we have used the measured values of the cosmological parameters [41] and assumed z, < 1 
in the last equality. The blazar emits 0(1 — 10 TeV) gamma rays into a jet or pair of jets [42]. In the 
simplest model of the jet, radiation is uniform within a cone of half opening angle and typically 
^jet — 5° corresponding to a solid angle of fijet — 0.024 sr. 

A gamma ray with energy ^ TeV at redshift z^^ is likely to scatter on optical and infrared 
extragalactic background light (EBL) and produce an electron-positron pair. The mean free path is 
given by where is the pair production cross section and driEBbi^^ ^ 77 )/de is 

the spectrum of EBL photons at redshift z^^. Assuming nEBL(^ 77 ) = (1 + ^ 77 )~^^ebl(^ = 0) it was 
argued in Ref. [40] that the mean free path can be reliably approximated as 

( 2 ^ 2 ) 

Ref. [40] estimates a range of values 0.3 < < 3 for the dimensionless coefficient, and we will take 

= 1 hereafter. The comoving mean free path is given by = (1 + As long as d^^ <C dg 

we can approximate and write the comoving mean free path as 
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After pair production the charged leptons acquire an energy ~ ^ 70 / 2 , and, to accuracy 
mejEe ^ 10“^ (the inverse boost factor), they travel in the same direction as the initial TeV gamma 
ray. The leptons produce secondary 7 -rays through inverse Compton (IC) scattering on cosmic mi¬ 
crowave background (CMB) photons. If Ecmb is the energy of a CMB photon at redshift then 
energy conservation gives the energy of the corresponding IC photon to be 
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where rrie — 0.511 MeV/c^ is the electron mass. When this photon reaches Earth (z = 0) it will 
have been redshifted to an energy of E^ = (1 + If the spectrum of the TeV blazar is 

known, one can use the spectra of the EBL and CMB to calculate the spectrum of GeV gamma rays 
arriving at Earth. Eor our purposes only the average relationships are required. At redshift the 
average energy of a CMB photon is (E^cmb (^77)) ~ ^CMB — (6 X 10 ^ eV)(l + 2 :^ 7 ), and the average 
energy of a GeV photon arriving at Earth is then 


Ej — ^ (1 -h Zjj 
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Evidently an initial spectrum of gamma rays from E^^ ^1 — 10 TeV are transferred into cascade 
photons with E^ 1 — 100 GeV energies. The scattering of leptons on CMB photons is a stochastic 
process that occurs with a typical mean free path /^fp. On average the leptons lose their energy via 
IC within the electron cooling distance 
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where (Jt — 6 . 65 x 10 “^^ cm^ is the Thomson scattering cross section and f/cMB — (l+2;^7)^(0.26eV/cm^) 

is the CMB energy density at redshift z^^. 

The halo emerges because TeV gamma rays directed off of the line of sight can induce a cascade 
that is deflected back toward the line of sight by the magnetic field. At the point of pair production 
(redshift z^^) let B and v be the magnetic field and the velocity of the lepton, respectively. Assuming 
that the coherence length of the magnetic field A |B|/| VB| is much larger than the electron or 
positron will probe an effectively homogeneous magnetic field. In this background the lepton follows 
a helical trajectory with gyroradius (Larmor radius) 

= where i? =, (2.7) 

c e B 


and — e is the charge of the electron. The component of v perpendicular to B is = v — (v • B)B 
where unit vectors are denoted by a hat. Using this expression we can write 


Rl = R ^1-(v-B)2 . 


( 2 . 8 ) 


If the magnetic field is frozen into the plasma, then its energy density redshifts like radiation, and we 
have |B| = with Bq the field strength today. Then we can we esti m ate 
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where we have also used 2 : 5 . 

By approximating throughout our analysis we assume that d^^/ds <C 1. Using 

Eqs. (2.1) and (2.3) this ratio is 
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where dg and are related by Eq. (2.1). We will primarily be interested in > 10 GeV and 
(is > 0.7 Gpc where the approximation is well-justified. Eor lower energy gamma rays or closer 
sources, the approximation begins to break down as the Earth is located inside of the developing 
cascade. 

It is also useful to compare the electron cooling distance and gyroradius R. Using Eqs. (2.6) 
and (2.9) and approximating 2 ;^^ 2:5 this ratio is 
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Since 27vR l is the circumference of the lepton’s orbit, this ratio indicates whether the lepton travels 
around the orbit many times De/R ^ 1 or whether it makes only a small arc 79e/i? <C 1. In the 
former case, called the “pair halo regime,” the cascade gamma rays are spread out over large angles; in 
the latter case, called the “magnetically broadened cascade regime,” the cascade photons are slightly 
spread out around the source. 


3 Morphology of the Cascade Halo 

We are interested in the size and shape of the cascade halo as it appears from Earth. In this section 
we first establish an analytic formalism for calculating halo maps, i.e. the orientation h(£^^) of GeV 
gamma rays reaching Earth. Next we introduce parameters that quantify the halo size and shape, 
which can be extracted from the halo map. 

3.1 Formalism and Assumptions 

We move between spherical, cylindrical, and Cartesian coordinates. The origin is located at the Earth 
and z is oriented along the line of sight to the blazar. The polar and azimuthal angles are denoted 
by (0, (/)), and the spherical and cylindrical radial coordinates are r and p. The sets of unit vectors 
{p, and {x,y,z} form right-handed orthonormal coordinate systems. It will be 

useful to note the relationships 

p = cos(/)x-h sine/)y x = cos(/)p - sine/) 

(p — — sin (/) X -h cos (py ’ y = sin p p -\- cos p p 

Eor plotting gamma ray arrival directions, it is convenient to introduce the lateral and transverse 
angular extent as 


^lat^Ocosp and T^trans = 6 ^ sin (/>. (3.2) 

The distinction between the lateral and transverse directions is arbitrary. Although other mappings 
are possible, this one has the convenient feature that the Euclidean norm is equal to the polar angle 

^ ^ '^trans- ^ gamma ray arriving at Earth is specified by an energy E^ and an orientation, 

which can be expressed as the pair ( 0 , p) or equivalently (^lat? ^trans) or fi = r. 
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We seek to study the development of the cascade semi-analytically by deriving a set of trigono¬ 
metric equations that can be solved (analytically if possible, numerically if not) to find the orientation 
of GeV gamma rays reaching Earth. In this sense our approach differs from a purely numerical sim¬ 
ulation in which TeV gamma rays are ejected from the blazar in random directions, the paths of the 
lepton and GeV gamma rays are calculated, and trajectories that do not intersect the Earth are dis¬ 
carded. As we will see, inspection of Eig. 1 leads to a set of three equations: the first arises from the 
trigonometry of the triangle, the second arises from the geometry of the lepton’s orbit, and the third 
ensures that the GeV gamma ray intersects with the line of sight. We solve these three equations for 
a given magnetic field configuration and gamma ray energy to obtain the orientation of the GeV 
gamma ray at Earth (0, 0) as well as the bending angle 5. This semi-analytic approach is applicable 
thanks to the following two well-justified approximations. 

Eirst, we assume that the lepton samples a homogeneous magnetic field, and thus its path is a 
simple helix. If A is the coherence length of the magnetic field, this condition is expressed as <C A. 
We saw in Eq. (2.6) that typically ^ 100 kpc, and since we will be interested in A > 1 Mpc this 
assumption is well-justified. If we were interested in smaller coherence scales, where the lepton 
motion is diffusive, then our approach would not be applicable. 

Second, we assume that the displacement of the lepton can be neglected. This ensures that, to 
a good approximation, the two gamma rays lie in the same plane and that they form the legs of a 
triangle as in Eig. 1. To verify that the longitudinal displacement along the field line is negligible, we 
need |v|| |T)e/|v| ds^d^Q\io ensure that the transverse displacement around the orbit is negligible 
we need Mm[De^ Rl] "C ds^ d^^. In Eqs. (2.1), (2.3), and (2.9) we saw that typically dg ^ 1 Gpc, 
^ 100 Mpc, and ^ 100 kpc. Then over all of the relevant parameters space we have 

De ^ 70 ’ ^he assumption is very well-justified. 

Although our approach only requires the above two approximations, we also use the following 
two assumptions as a matter of convenience. Eirst, it is important to remark that many aspects of 
the cascade are stochastic in nature: the spectrum of TeV gamma rays emitted by the blazar, the 
distance traveled by the TeV gamma rays before pair production, the spectrum of EBL photons, the 
distance traveled by the leptons before IC, the number of IC photons emitted before electron cooling 
becomes appreciable, and the spectrum of CMB photons that are up-scattered. Eor simplicity we 
neglect the stochastic spread in each of these various parameters, and we fix them equal to their 
average values given in Sec. 2. As a result, we obtain a deterministic relationship between the energy 
of gamma rays reaching Earth and their orientations on the sky; we call this function the halo map 
These halo maps are useful tools for studying the connection between halo morphology and 
the underlying magnetic field since they can be calculated quickly, without sophisticated numerical 
simulation, and they capture the characteristic features of the halo size and shape. It is expected that 
properly taking account of the stochasticity will lead to a significantly “smeared” version of the halo 
maps shown here, since the variance in the random parameters is typically 0(1). We discuss the 
stochastic smearing further in the conclusions. Sec. 6, as a direction for future work. 

As a second simplification, we neglect the blazar’s jet structure and assume that emission from 
the blazar is isotropic. If we were to properly treat the angular distribution of radiation from the blazar, 
then only a portion of our halo maps would be visible. We return to this point in the conclusions. 
Sec. 6, where we also show a few halo maps that have been calculated with the jet restriction in place. 

3.2 Constraint Equations 

We obtain the first in the set of three equations by applying the law of sines to the triangle in Eig. 1: 

sinO = sin5 . (3.3) 

ds 
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Here 9 is the polar angle that a GeV gamma ray arriving at Earth makes with the line of sight to the 
blazar, and the bending angle 6 is the angle between the orientation of the initial TeV gamma ray and 
the final GeV gamma ray. The second equation is a relationship for the bending angle 6: 

1 - cos (5 = (^1 - (vi • B)2) (l - cos(De/i?)) . (3.4) 

To derive Eq. (3.4), we write the magnetic field at the point of pair production as B = S n||, and we 

write the initial lepton velocity as ny + n± where fiy • = 0. Then ^’y = • B and 

^y c^. The lepton velocity at time t after pair production is 

v(f) = ?;y ny + v± cos{(jjt) fi^ =F v± ny x (3.5) 

where u = vj_/Rl = c/i? is the angular frequency of the orbital motion. The sign ambiguity in 

the last term is related to the charge of the lepton; the (—) is for positrons and (+) is for electrons. 
The lepton travels an electron cooling distance^, in time r = De/c. Since the gamma rays 
are approximately tangential to the lepton trajectory, the orientation of the initial TeV gamma ray 
is v(0), the orientation of the final GeV gamma ray is v{De/c), and the bending angle satisfies 
cos 6 = v(0) • v{De/c), which gives Eq. (3.4). 

Before discussing the third equation, it is useful to consider the limits of small and large lepton 
deflection. Recall that D^/R — ccr is the angular deflection of the lepton as it travels a distance 
around the gyro-circle, and we saw in Eq. (2.11) that D^jR (x Eor sufficiently weak 

magnetic field or high gamma ray energy we have D^/R <C 1 and Eqs. (3.3) and (3.4) reduce to 

e - (Vj • B)2 0ext and 5 « - (v* • B)2 E ( 3 . 6 ) 


where 


Oex. = . (3.7) 

dsR 

In ©ext we find the familiar expression for the angular extent of the halo [14]. Using the numerical 
estimates from Sec. 2 we have 


©ext ^ (0.68°) 


Bo 


10-15 G 


lOGeVy VlGpc) V 1-24 ) 


(3.8) 


Since ©ext the angular extent of the halo grows larger as the field strength is increased, which 
is characteristic of the MBC regime. 

In the opposite regime where the lepton deflection is large we have D^/R ^ 1, and the lepton 
makes multiple orbits before IC occurs. If (v^ B )2< 1/2 then Eq. (3.4) has a discrete set of solutions 

in) 

at energies ^ where the halo reaches a maximum angular extent 

0 «©max and ^ ~ I (3.9) 

^Throughout our analysis we assume that the lepton always travels an electron cooling distance, and that IC scattering 
only occurs once, at this point. In reality the mean free path of the lepton is shorter than De and multiple IC photons 
are emitted. These multiple emissions are an example of the stochastic effects, discussed previously, that we neglect for 
simplicity. We perform a preliminary investigation of the stochasticity in Sec. 6.” 
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with 


©max = arcsin 



~ (10°) 


4”^ V 4 + "^ 

10 GeV j V 1-24 
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In this limit, known as the PH regime, the angular extent of the halo is not proportional to the magnetic 
field strength, but instead it is restricted only by the geometry of the TeV gamma ray propagation. 

We finally turn to the third constraint equation, which needs to enforce that the TeV and GeV 
gamma rays approximately lie in a plane of constant 0, as shown in Fig. 1. To ensure that the charged 
lepton is not deflected out of this plane, the Lorentz force F = (e/c) v x B must be normal to 0. 
Of course, as the lepton follows its helical path, the direction of its velocity changes and so too does 
the direction of the Lorentz force. Then we must average the Lorentz force over the trajectory of the 
lepton. 

Let x(f) be the helical trajectory of the lepton, v(f) = rfx/rff be its velocity, be the time 
of pair production, and r = D^/c be the time elapsed before IC up-scattering. Neglecting the 
cosmological expansion, which is not relevant on such short time scales, the impulse imparted on 
the charged lepton at redshift is given by 


J = 




df v(f) X B(x(f), f) 


(3.11) 


where the ± is related to the charge on the lepton. If the magnetic field is static and homogeneous 
over the path of the lepton, we can pull it out of the integral, and the integrand contains only v(f). 
Defining the time averaged electron velocity as 


1 

Vavg^-y dt\{t) , 


(3.12) 


we can write the impulse as 


J = i- Vavg X B . (3.13) 

Since the magnetic field does no work, the magnitude |v(f)| c is fixed. Then using the geometry 

shown in Fig. 1, we see that Vavg must bisect the angle 5, and it can be written as 

-6 >)z. (3.14) 


-9 


p — cos 


Writing also 


B = hpp + + bzZ 


(3.15) 


we have 


^ avg 


X B • ^ = —bp cos ( f 0 


bz sin { 2 “ I 


(3.16) 


Typically bp or bz will depend on the azimuthal angle cj), and then this equation fixes the plane (normal 
to 0) in which lie the two gamma rays and the line of sight. 
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Now we summarize the three constraint equations. To ensure that the motion remains in the 
plane, we impose Eq. (3.16): 

bp cos sin - 6*^ = 0 (3.17) 

where Eq. (3.15) gives the decomposition of the magnetic field into cylindrical coordinates. Erom the 
geometry of the gamma ray trajectories, we have the law of sines in Eq. (3.3), 

sinO = sin5 , (3.18) 

dg 

and finally the bending angle is given by Eq. (3.4), 

1 — cos 5 = ^1 — (bpsm{6 — 9) — bzCos{6 — 9))^^ ^1 — cos(De/di)^ , (3.19) 

where we have written the initial lepton velocity as 

= sm{S — 9)p — cos(5 — 0)z . (3.20) 

The constraints in Eqs. (3.17)-(3.19) can also be derived from the single vector equation 

d^^Vi + Ax + Lvf + dgZ = 0 , (3.21) 

which ensures that the cascade photon reaches Earth. Here L is the distance from the IC scattering 
to the observation point and Ax is the displacement of the lepton between pair production and IC 
scattering. To find Ax, we can integrate Eq. (3.5). However, in the limit that <C dg, the Ax 
term can be dropped from the equation. Trigonometric manipulation of Eq. (3.21) then once again 
leads to the above constraint equations. 

3.3 Shape Parameter 

In the next section we consider various static magnetic field configurations B(x), and we solve 
Eqs. (3.17)-(3.19) for 0, 0, and 6 to determine the halo map n(£’^). Having solved Eqs. (3.17)-(3.19) 
we construct the halo map as the radial unit vector 

n{Ery) = sin 9 cos 0 x + sin 9 sin 0 y + cos 9 z (3.22) 

that points from the Earth toward the arriving gamma ray of energy E^. Note that n{Ej) may be 
multi-valued meaning that gamma rays of a particular energy may appear from multiple directions in 
the sky. All of the information about the halo size and shape is contained in the function h{Ej). Eor 
instance, the size of the halo is given by &{E^) = arccos[n(£’^) • z]. 

There are many ways to quantify the halo shape and orientation. Since we are interested in 
probing magnetic helicity, we are motivated to consider a parity-odd Q-statistic [36, 37], which is 
sensitive to the sign of the magnetic helicity. The statistic is defined as the triple product of vectors 
sampled from the halo map n{Ej) at three different energies. Eor illustrative purposes, we consider 
a few energy combinations: 


Qio, 30,50 = nio X n3o • nso ( 3 . 23 a) 

Q25,30,35 = ^25 X 1130 * 635 ( 3 . 23 b) 

Q8,io, 12 = ns X fiio • ni2 ( 3 . 23 c) 

Q38,40,42 = n 38 X h^Q • 1142 ( 3 . 23 d) 
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where fi^ is shorthand notation for n{E^ = # GeV). We will see that different energy combinations 
are sensitive to different ranges of parameters, namely magnetic field strength and coherence length. 
Since may be multi-valued in general, one can calculate Q by first averaging over the multiple 

arrival directions for a given energy. 

It is also useful to write the triple product Qabc = x • Ac as 

Qabc = sin dab sin ipabc (3.24) 

where 0 < dab < tt is the angle between and n^, and — 7 r /2 < Lpabc < 7 r /2 is the angle between 
Ac and its projection onto the plane spanned by A^ and A^. If A^ and A^ are collinear then dab = 0 
and Qabc vanishes; if Ac lies in the same plane as A^ and A 5 then cpabc = 0 and Qabc also vanishes. 
The sign of Qabc is controlled by the sign of (pabc^ which depends on whether Ac is “in front of” or 
“behind” the plane normal to A^ x A5. Since it is odd under reflections one sees that Qabc is a measure 
of parity violation in the halo map. 

4 Halo Morphology for Specific Magnetic Field Configurations 

In this section we study the size and shape of the GeV halo for various specific magnetic field configu¬ 
rations. The configurations we consider are simplified and do not realistically model the intergalactic 
magnetic field. However, these examples serve to illustrate the parametric relationships between the 
field configuration and the halo morphology. Previous studies have focused on the size information 
alone, and we will see that the shape information provides insight into the magnetic field’s orienta¬ 
tion and helicity. We consider five different non-helical and helical magnetic field configurations, as 
shown in Fig. 2. 

4.1 Case 1: Uniform Magnetic Field Parallel to Line of Sight 

The simplest configuration is a homogeneous magnetic field oriented along the line of sight with the 
blazar, 

B = -z . (4.1) 

The three constraint equations, Eqs. (3.17)-(3.19) , reduce to 

sm{6/2 - 9) =0 (4.2a) 

sin (9 =^sm6 (4.2b) 

ds 

1 — cos 5 = sin^(5 — 0)^1 — cos{De/R)^ . (4.2c) 

There is a trivial solution with 6 = 9 = 0 corresponding to gamma rays oriented along the line of 
sight, V = —z, that are not deflected by the magnetic field. There is also a nontrivial solution, 

s = 20 = 2 cos~^ ^ E = (^2n+l)7r (4.3) 

where n is an integer. Since the azimuthal angle (j) does not appear in these equations, the solution will 
be rotationally synunetric about the line of sight to the blazar. Also, since Da/R and d^^/dg depend 
on energy, there will only be a discrete set of energies for which a solution exists. The solution can 

be understood in physical terms: the velocity component along the magnetic field remains constant, 

and the velocity component perpendicular to the magnetic field gets reflected, and so the triangle in 
Fig. 1 is an isosceles triangle. Apart from these solutions, the magnetic field deflects other gamma 
rays away from the line of sight, and they do not reach Earth. 
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CASE 1 


CASE 4 



Figure 2. The five magnetic field configurations that we consider. Cases 1-3 are homogeneous field configura¬ 
tions that have different orientations with respect to the line of sight with the blazar. Cases 4-5 are helical field 
configurations with their wavevectors oriented either along or normal to the line of sight. 


4.2 Case 2: Uniform Magnetic Field Normal to Line of Sight 


Next we consider a homogeneous magnetic field that is oriented normal to the line of sight with the 
blazar (see Fig. 2). Without loss of generality we can align the Cartesian coordinate system with the 
magnetic field such that 


B = y = sin (/) p + cos (j) cj) . 

(4.4) 

and Eqs. (3.17)-(3.19) reduce to 



sin(/) cos(5/2 — 0) = 0 


(4.5a) 

sin 9 = sin 6 

dg 


(4.5b) 

1 — cos 5=^1 — sin^(5 — 9) sin^ 

4>J (^1 - cos{De/R)J . 

(4.5c) 

For a given ganuna ray energy there is a solution 



(/) = 0,7r , sin0=^^sin5 

dg 

, and cos 6 = cos 

R 

(4.6) 


where enters through d^^/dg and De/R, see Eqs. (2.10) and (2.11). Recall that De/R > 0 is 
unbounded from above but 0 < 5 < tt and 0 < 0 < 7r/2. For this magnetic field configuration, the 
trajectories of all the gamma rays lie in the 7/ = 0 plane where 0 = 0, tt. In the limit of small lepton 
deflection, DelR<^ 1, the solution further simplifies to 


e 




©ext — 


d^^D^ 

dgR 


(4.7) 
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as in Eqs. (3.6) and (3.7). 

The halo map corresponding to the solution in Eq. (4.6) is shown in Eig. 3. We have mapped 
the polar and azimuthal angles (0, (/)) into the lateral and transverse angular extent (^lat? ^trans) using 
Eq. (3.2). The ^-component of velocity is conserved for B = y, and only gamma rays lying in 
the plane y = 0 normal to the magnetic field are deflected back toward the line of sight. The halo 
map is double-valued as there are precisely two gamma ray arrival directions for each energy. 

Along one branch the electron was deflected back toward Earth and generated the secondary IC 
gamma ray, and on the other branch it was the positron. 

The energy dependence of the polar angle is shown in Eig. 4 for various parameter combinations. 
The highest energy gamma rays experience the smallest deflection, and they arrive closest to the line 

of sight (9 = 0). In this regime, the small angle approximation is valid, and we have the scaling 

— 3/2 

9 ^ ©ext from Eq. (3.7). Lower energy gamma rays are found farther from the line 

of sight, and there is an energy gradient. As the energy decreases further, the bending angle 8 = 
De/R oc l/E^ continues to grow, as per Eq. (2.11), and eventually 6 = tv 12 where the electron 
experiences a deflection of 90°. This corresponds to an energy 

” - (»■« G'V) (is^) (^)". <‘''® 

which also serves to indicate where the small bending approximation breaks down. Since sin 5 = 1 
is maximized at this point, the halo achieves a maximum angular extent (cf. Eq. (3.10)) 

e„„ = atcsin = (59«) ' • (4.9) 

Still lower energy gamma rays arrive closer to the line of sight, because the lepton is bent more than 
90°. 

Eor this case the halo map is wide in lateral extent and narrow in transverse extent as seen 
in Eig. 3. In principle one could measure the orientation of the magnetic field by measuring the 
orientation of the halo map. This is an example of how shape information can probe additional 
aspects of the IGME beyond just its field strength. Of course, we have assumed that the magnetic 
field is uniform over the scale probed by the TeV gamma rays, Le. if the magnetic field coherence 
length is A, then we have implicitly assumed X ^ ^ 100 Mpc. In a realistic setting it is more 

likely that the magnetic field forms domains smaller than . If the magnetic field is statistically 
isotropic across domains, then different leptons probe random orientations of the field, and the halo 
will resemble a more familiar, rotationally symmetric halo map. Apart from this isotropization, the 
discussion of this section is largely unchanged, and specifically Eq. (4.6) still gives the relationship 
between the gamma ray energy and polar angle. 

4.3 Case 3: Uniform Magnetic Field with Arbitrary Orientation 

Next we consider a homogeneous magnetic field that has a component along the line of sight to the 
blazar. This is a generalization of the previous two cases. Without loss of generality we can write the 
magnetic field configuration as 

B = (cos /3 y — sin /3 z) = cos /3 sin 0 p + cos /3 cos 0 ^ — sin /3 z (4.10) 
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Figure 3. The halo map for Case 2. Gamma rays arrive at Earth collimated into a line with zero transverse 
extent. The black dots indicate = 100, 20,15,10, and 5 GeV with the higher energy gamma rays arriving 
near the origin. The gray dashed lines are curves of constant polar angle 0 . 






Figure 4. Angular extent of the halo as a function of gamma ray energy for Case 2. Each panel shows the 
exact solution (black line) from Eq. (4.6), the small bending approximation (red dashed) from Eq. (4.7), and 
the geometric limit (black dashed) from Eq. (3.10). The top-right panel shows the same parameters as in Eig. 3. 
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Figure 5. The halo map for Case 3. If we had taken < 0 then the halo maps would be reflected in the vertical 
direction. The black dots and gray dashes have the same meaning as in Fig. 4. 


where the skew angle /3 controls the component of B along the line of sight. Eqs. (3.17)-(3.19) reduce 
to 


siiKp = tan/3 tan((3/2 — 9) (4.11a) 

sin 9 = sin S (4.1 lb) 

dg 

l_c„s{= (l-sm2,3-J?i|§^)(l-cos(D,/iS)) (4.11c) 

where we have used the equation for sin (j) to simplify the third equation. 

Numerically solving Eq. (4.11) leads to the halo maps shown in Eig. 5. In comparing with 
Case 2 from Sec. 4.2 we see that the halo map is no longer restricted to a line, but instead the halo 
acquires a transverse extent. We vary the skew angle /3 in the left panel of Eig. 5. In the limit that 
/3 goes to zero, we regain the line-like halo map of Case 2, and in the limit that (3 goes to 90°, we 
regain point-like halo map of Case 1. We vary the field strength in the right panel of Eig. 5. In the 
PH regime where Sq < 10 x 10“^^ G, the halo size is proportional to the field strength (green and 
magenta curves), while in the MBC regime where Bq > 10 x 10“^^ G, the halo size is limited by 
the geometry (brown and teal curves). 

Our analytic approach to calculating the halo map has an advantage over numerical shooting 
techniques insofar as we can solve for the halo map analytically in certain limiting regimes. To 
demonstrate this point, first consider the limit of small skew angle /3 <C 1 in which Case 3 reduces to 
Case 2. In this limit, Eq. (4.11) becomes 


(j) ^ (3 tan(5/2 — 9) 


9 ~ arcsin 


^70 


Sin 


R 


and 


5 




De 

R ’ 


(4.12) 
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and it follows that 


f^iat « arcsin ^ ^ ^ 

■t^trans ^ ^ arcsin ^ arcsin 


This behavior is seen in the left panel of Fig. 5. Second, consider the small bending angle regime 
De/R <C 1 (weak magnetic field) where Eq. (4.11) becomes 


sin (3 Dg 


e 


G. 


ext 


and 6 ^ cos /3 , 
R 


2 cos^ /3 i? ’ ^ cos [3 

and ©ext Bq was given by Eq. (3.7). Then the lateral and transverse extents are 


'^lat ~ 
'^trans 


©ext 
COS /3 

^ sin /3 Dq ©ext 
2 cos^ /3 R cos /3 


(4.15) 


(4.16) 

(4.17) 


Observe that ^trans is suppressed with respect to 7?iat by an additional factor of {De/R) oc Bq. 

We quantify the halo shape using the triple product Q-statistic, given by Eq. (3.23). Since the 
magnetic field configuration under consideration is not helical, we expect Q = 0. In fact this is 
immediately evident from the symmetry of Eig. 5: the gamma rays on the branch in the first quadrant 
contribute Q < 0 while those in the second quadrant contribute Q > 0, and upon summing the 
two branches, they cancel. However, this cancellation is possible in part because we have assumed 
isotropic emission from the blazar. In practice, the jet may only illuminate a small patch of the halo 
map. Thus, to demonstrate the parametric dependence and typical scale of Q it is illustrative to 
calculate the statistic using only the gamma rays in one of the two branches. 

We evaluate (5io,30,50 from Eq. (3.23) and show the results in Eig. 6. Similar results are obtained 
for the other energy combinations, and we do not show them here. The parameters are chosen to 
correspond with the halo maps in Eig. 5. The statistic Qio,30,50 becomes small (i) in the limit (3 
0 where the halo map approaches a straight line, (ii) in the limit (3 90° where the halo map 

approaches a point, and (hi) in the limit So ^ 0 where the angular extent of the halo decreases. 
Increasing the field strength grows (5io,30,50 until the crossover from the MBC to the PH regime at 
So ^ 50 X 10“^^ G. Eor larger So the statistic first decreases and then begins to oscillate, similar to 
the behavior of the halo size, seen in Eig. 4. 


4.4 Case 4: Helical Magnetic Field with Wave Vector Parallel to Line of Sight 

We now turn our attention to helical magnetic field configurations. The simplest configuration con¬ 
sists of a single circular polarization mode with wavelength A and wavevector k = (27r/A)z oriented 
along the line of sight with the blazar: 

B = cos('0 + 2'kz/\) y + (7 sin('0 + 2'kz/\) x . (4.18) 

The spatial coordinate z should not be confused with the redshift The three parameters are the 
coherence length A, the handedness index cr = ±1, and the phase shift The handedness index 
controls the sign of the magnetic helicity density, 

H = B-VxB = a^|Bp. (4.19) 

A 
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Figure 6. The magnitude of the parity statistic, given by Eq. (3.23), for Case 3. If we had averaged over all 
gamma rays, we would find Qio, 3 o ,50 = 0 since the gamma rays in the first and second quadrants of Fig. 5 
cancel. Here we just show Qio, 3 o ,50 for the gamma rays in the second quadrant. 


The case a = 0 corresponds to a non-helical, linearly polarized plane wave, while a = +1 corre¬ 
sponds to left-circular polarization and —1 to right. In the subsequent analysis, one should bear in 
mind that varying A at fixed |B| implies that the magnetic helicity is being varied. 

It is convenient to define the angle /3eff = ^ + 2tizI\. Then using the geometry of Fig. 1 the 
longitudinal coordinate at the point of pair production is z = cos(5 — 9), and 

Ptt 

^eff (5, 6) ='ll) + — {ds - djg cos{5 - 9)) (4.20) 

is the effective skew angle. For this case, Eqs. (3.I7)-(3.I9) reduce to the set of equations 


tanc/) = —a 

• n ^70 

Sin 9 = —^ 

dg 

cos 6 = cos 


tan/3eff((5,6') 
sin 5 

R ’ 


(4.21a) 

(4.21b) 

(4.21c) 


which can be solved analytically. The solutions are shown in Fig. 7. The most striking feature in 
these figures is that the halo map forms a spiral pattern. The handedness of the spiral is controlled by 
the helicity of the magnetic field, parametrized here by a = ±1. As the phase shift ^ is varied, the 
halo map is uniformly rotated clockwise or counterclockwise. 

As we vary the coherence length A the spiral becomes flatter or tighter. In the limit of large 
coherence length X d^^ ^ 100 Mpc, the cascade takes place in an effectively homogeneous 
magnetic field, B ^ y. Then we regain the behavior of Case 2 from Sec. 4.2 in which the gamma 
rays propagate in a plane and arrive at Earth collimated into a line with small transverse extent. In 
the opposite limit of small coherence length, A <C 100 Mpc, the TeV gamma rays sample the 

magnetic field at a random phase. In terms of the halo map, this translates into a tightly wound spiral 
with multiple cycles. Varying the magnetic field strength has the same effect as in the previous cases. 
In the PH regime where the field is weak, the angular size of the halo grows with increasing field 
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Figure 7. The halo map for Case 4. We have taken a = +1 and ^ = 0 . For a = —1 the handedness of the 
spiral is reversed, and for a phase shift 7 ^ 0 the halo map is uniformly rotated, but otherwise the structure 
remains unchanged. In the left panel, if A is further decreased below 100 Mpc the spiral becomes tighter, and 
if A is further increased above 1 Gpc the halo map asymptotes to a straight line. 


strength, while in the MBC regime where the field is strong, the angular size of the halo is limited by 
the geometry (brown curve in right panel). 

The halo maps in Fig. 7 display a clear parity-violation that should be captured by the Q- 
statistics. However, the halo map is double-valued, Le. there are two gamma rays at each 

energy, which correspond to the two branches of the spiral in Fig. 7. If we first average n{Ej) 
over its multiple solutions, we would obtain Avg[n(£^^)] = 0 and therefore Q = 0. This result is 
a consequence of the high symmetry of the system under consideration: we have assumed that the 
wavevector of the magnetic field is oriented along the line of sight, and we have allowed for isotropic 
emission from the blazar. As we discussed for Case 3 of Sec. 4.3, under more realistic conditions the 
blazar’s jet will only illuminate a part of the full halo map, and then the cancellation is disrupted. To 
model this effect, we calculate Q along a single branch of the spiral. 

We calculate the Q-statistics using Eq. (3.23) and show the parametric dependence on coherence 
length A and field strength Bq in Fig. 8. The statistic has the same qualitative behavior for each energy 
combination: (i) |Q| decreases for small Bq, (ii) \Q\ decreases for large A, and (hi) Q oscillates 
rapidly from positive to negative values for small A. When the magnetic field is weak, the halo has 
a small angular extent, as shown in the right panel of Fig. 7, and Qabc = x • he decreases as 
ha, h^, and he become approximately collinear (see also Eq. (3.24)). When the coherence length is 
large the halo map resembles a straight line, as seen in the left panel of Fig. 7, and Qabc decreases as 
ha, h^, and he become approximately coplanar. In this regime, the magnetic field appears uniform, 
and therefore non-helical, on the scale probed by the TeV gamma rays. When the coherence length 
is small the halo map resembles a tight spiral, as seen in the left panel of Fig. 7, and the three gamma 
rays used to construct Q may not lie on the same cycle of the spiral. In this case, Q may take positive 
or negative values depending on the energies at which the spiral is sampled, and as A decreases further 
and the spiral becomes more tightly wound, the sign of Q oscillates. We will use Aqsc to denote the 
coherence length at which Q begins to oscillate. 
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Figure 8. The parity statistics, given by Eq. (3.23), for Case 4. For large A the sign of the Q-statistic is given 
by sign[(5] = —a = — sign[H]. As the coherence length is lowered, the halo map begins to spiral around 
the line of sight, as seen in Fig. 7, and the sign of Q oscillates. Statistics calculated from different energy 
combinations have the same qualitative behavior, but they differ in the magnitude of Q and the value of A at 
which the oscillations begin. 


The four cases in Fig. 8 display important quantitative differences between Q-statistics con¬ 
structed from different energy combinations. Comparing the bottom two panels, IQs, 10,12! > |Q38,40,42 
we see that lower energy gamma rays leads to a larger value for |Q|. This is simply because lower 
energy gamma rays are more easily deflected and lead to a larger halo in the MBC regime, see 
Eq. (3.8). Comparing the top two panels, Aosc({10, 30, 50} GeV) > Aosc({25,30, 35} GeV), we 
see how the energy spacing affects the coherence length scale below which Q begins to oscillate. 
As A decreases and the halo map becomes a more tightly wound spiral, and the gamma rays from 
more closely spaced energies, in this case {25, 30,35} GeV, remain on the same cycle of the spiral 
longer than more widely spaced energies. In the two lower panels the energy spacings are identical 
but nevertheless Aosc({38,40,42} GeV) < AoscdS, 10,12} GeV). This is because the spiraling 
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behavior becomes more pronounced as the gamma ray energy is lowered; see Fig. 7 and note that the 
separation between the 5 and 10 GeV gamma rays is much larger than the separation between the 10 
and 15 GeV gamma rays. 

4.5 Case 5: Helical Magnetic Field with Wave Vector Normal to Line of Sight 

Finally we consider a more generic helical magnetic field configuration. We do not require that the 
wave vector is oriented along the line of sight, but for simplicity we don’t allow it to have a general 
orientation either, and instead we require that k is normal to the line of sight with the blazar. Without 
further loss of generality we can write k = (27r/A)x. The magnetic field configuration is 

B = cos('0 + 27rx/A) y — (J sm{^ + 2'kx/\) z 

= cos ('0 + 271x1 X) sin 0/0 + 008(0 + 27rx/A) cos 00 — a sin (0 + 27rx/A) z . (4.22) 

In this case the phase shift 0 plays a nontrivial role. For specific values 0 = 0 , 7 r/ 2 , 7 r,-- - the 
magnetic field is either symmetric or antisymmetric when reflected across the line of sight {x —x 
and 0^—0), but for general 0 there is no such symmetry. Once again the handedness index a 
controls the magnetic helicity via Eq. (4.19). 

It is convenient to introduce the effective skew angle 

0eff(0 0,0) = 0 + sin(5 — 6) cos 0 (4.23) 

A 

where we have used x = sin(5 — 9) cos 0. Then Eqs. (3.17)-(3.19) reduce to 

sin0 = (j tan 0 efr (^7 6^5 0 ) tan (5/2 — 9) (4.24a) 

sin 9 — sin 6 (4.24b) 

dg 

1 - cos 5 = (l - sin^ ^eff 0,4>) e ) ) ~ cos(L>e/^)) • (4.24c) 

These equations have a rich and interesting family of solutions, but as a result they cannot be solved 
analytically as in Case 4. Instead we solve Eq. (4.24) numerically. 

The halo maps are shown in Fig. 9, and we have chosen the same parameters as in Fig. 7 
to facilitate comparison with Case 4. In each of the three panels we have taken the handedness 
index a = +1, and the associated parity-violation is evident in the “S”-like shape of the halo maps. 
Choosing a = —1 reflects the halo maps across the vertical axis, which flips the handedness. 

In the first panel, the coherence length is reduced from A = 1000 Mpc to 50 Mpc. The behavior 
in the large A regime is the same as in Case 4: when the coherence length is much larger than the 
scale of the cascade, d^^ ^ 100 Mpc, the gamma rays probe an effectively homogeneous magnetic 
field, and we regain the line-like halo map that was originally seen in Case 2 of Sec. 4.2. The small A 
behavior is distinctly different than in Case 4 where we encountered a spiral-shaped halo map. Now 
the wavevector crosses the line of sight to the blazar. As A is reduced, the surfaces of constant phase 
become compressed in the lateral direction, and the “S”-like halo map becomes “squeezed.” 

In the second panel, the magnetic field strength is increased from Sq = 2 to 16 x 10“^^ G. In 
the small Bq regime, we once again regain the behavior of the previous cases: a smaller field strength 
translates into a smaller halo. In the large Bq regime, on the other hand, a new phenomenon emerges: 
the halo map acquires multiple disconnected branches. In this panel only, all of the black dots denote 
Ej = b GeV, and we see that for Bq = 8 and 16 x 10“^^ G the halo map n{Ej) has six distinct 
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Figure 9. The halo map in Case 5. We have taken a = +1 in all panels, and for a = —1 the halo maps are 
reflected across the vertical axis. In the first and second panels we fix t/; = 0. In the first and third panels, 
the black dots have the same meaning as in the previous halo maps: they indicate = 100, 20,15,10, and 
5 GeV. In the second panel, the black dots show only = b GeV, the black square indicates E^ = 7 GeV 
for Bq = 8 X 10“^^ G, and the black triangle indicates E^ = 11 GeV for Bq = 16 x 10“^^ G. 


values. In the previous cases, the halo map was only double-valued. The new branches contain only 
low energy gamma rays, which are more easily deflected, and the bifurcation points are denoted by a 
black square and black triangle for Bq = 8 and 16 x 10“^^ G, respectively. As the field strength is 
further increased or the coherence length lowered, additional branches will emerge. 

In the third panel, we vary the phase parameter ^ from 0 to tt. If ^ (0 mod tt) then the 

reflection symmetry of the halo map is disrupted. This was not the situation in Case 4 where varying 
ijj simply lead to a uniform rotation of the spiral-Iike halo map. 

To quantify the parity-violating features in the halo maps, we calculate the Q-statistics using 
Eq. (3.23) and show the results in Fig. 10. As we have seen, the halo map becomes multi-valued 
for large Bq and small A, and it displays multiple branches. To generate Fig. 10 we calculate Q 
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Figure 10. The parity statistics for Case 5. We calculate Q using only the main solution branch. The plots are 
truncated at small A and large Bq where numerical issues arise. For small A and large Bq where there can be 
multiple branches, this approach probably underestimates Q. 


using only the main branch, which is continuously connected to the line of sight. By neglecting the 
branches with larger angular separation from the line of sight, we would presumably underestimate 
Q in this regime. 

Over the parameter space shown Fig. 10 we have sign[(5] = —a = —sign[7{], where T-L > 0 
for the left-circular polarization mode. In the limit of large A we have the same behavior as in Case 
4, see Fig. 8. In the limit of small A the Q-statistics decrease because the main branch is squeezed, 
and the angular extent of the halo is smaller. Also in comparing with Case 4, we see that the scaling 
with Bq is different: in Fig. 8 we found \Q\ ^ Bq but Fig. 10 implies that \Q\ ^ Sq for Case 5. As 
such, the magnitude of Q very quickly decreases with decreasing field strength. 
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5 Implications for Helicity Measurement 

By considering particular realizations of the IGMF, our analysis reveals the parametric behavior of 
the halo size and shape on the various magnetic field variables, specifically the field strength Sq, 
the coherence length A, and the magnetic helicity density T-L = zb(27r/A)SQ. In the “corners” of 
the parameter space, where either or A is very large or very small, the halo map behaves in a 
way that may make a measurement of the magnetic helicity more challenging than in the “central” 
parameter regime. For instance, small So implies a small halo, which may not be distinguishable 
from a point source with a telescope’s finite angular resolution. In this section we will demarcate the 
various parametric regimes and discuss the challenges posed by each. 

The quantitative analysis in this section uses the simplified field configurations considered in 
Sec. 4, but we expect qualitatively similar results (parameter space boundaries) even if the interven¬ 
ing magnetic field configuration is not of one of the forms we have discussed, e.g. if the field is 
stochastically homogeneous and isotropic. The inclusion of stochastic variables in the development 
of the cascade, e.g. spectrum of the EBL or scattering probabilities, will smear out the halo patterns 
we have seen, and non-cascade photons introduce noise in our signal and dilute the halo. A rigorous 
evaluation of experimental sensitivities, even for a given set of experimental parameters, will require 
more information on blazar sources (gamma ray flux and spectrum, jet orientation and structure) and 
the background noise. 

Strong Field Regime 

For a strong magnetic field, at energies such that the gyroradius R is smaller than the typical distance 
traveled by the charged lepton see Eq. (2.11), we have D^/R > 1. This has two consequences. 
First, the maximum angular extent of the halo is no longer tied to the magnetic field strength, but 
instead it is fixed by the geometry as in Eq. (3.10). Then if a halo is seen, one can infer the presence 
of an IGMF, but one cannot measure the field strength from the halo size alone. Second, since the 
charged leptons can make a complete orbit around the gyrocircle, the GeV gamma rays will be emitted 
isotropically. This reduces the flux by roughly f]jet/47r, where fijet is the solid angle of the blazar’s 
jet, and makes it more difficult to see the halo. 

Therefore in the strong field regime, measurements of the cascade halo have a reduced capacity 
to probe the parameters of the IGMF. To avoid this regime, we require 



(5.1) 


which leads to an upper bound on the magnetic field strength. 



(5.2) 


upon using the expression for D^/R in Eq. (2.11). We plot this boundary in Fig. 11. It divides the 
PH regime D^IR> \ from the MBC regime D^/R < 1. 

Weak Field Regime 

A gamma ray telescope cannot pinpoint the arrival direction of a gamma ray with arbitrary precision. 
If the angular resolution is too poor, then the halo cannot be distinguished from a point source, and 
one cannot use halo size and shape measurements to probe the IGMF. This issue becomes especially 
relevant for weak magnetic fields, which do not induce much bending and lead to a smaller halo {cf. 
Eq. (3.8)). 

Angular resolution is quantified by the point-spread function (PSF), which is the probability 
distribution function for the angle between the true and reconstructed arrival direction of a gamma ray 
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of given energy. The 68 % confinement radius is the angular radius containing 68 % of the probability. 
The confinement radius of the Fermi-LAT (Large Area Telescope) [43] is well-approximated by the 
following empirical formula:^ 


69qs{E^) - ( 0 . 11 °) 



E, \ 

7.9 GeV ) 


- 1.62 


(5.3) 


Although our analysis is not specific to the Fermi-LAT instrument, we use this confinement radius as 
fiducial point of reference. 

In the weak field regime, the angular extent of the halo is given approximately hy 9 ^ ©ext (£’ 7 ) 
from Eq. (3.8), regardless of the specific field configuration under consideration. If ©ext is sufficiently 
large compared to 69qs, then the halo can be distinguished from a point source and its angular extent 
measured. In fact, if the detector response is known very well (negligible systematic error), then halos 
as small as 59q ^/can be probed when a large number of halo gamma rays are visible. Thus 
we assess when the telescope will be able to distinguish the halo from a point source using 


©ext(£^7) ^ 


89qs{E^) 


(5.4) 


where is the number of gamma rays collected at energy E^. 

The requirement of sufficient angular resolution in Eq. (5.4) leads to a lower bound on the 
magnetic field strength. 


^ (0.16 X 10-1^ G) 



E, \ 
7.9 GeV J 



E^ 

10 GeV 



/ I -\- Zs 

V 1.24 


3 

(5.5) 


which is shown in Fig. 11. For smaller Bq the cascade halo is too small to distinguish from a point 
source given an angular resolution comparable to the Fermi-LAT We show a scenario with small 
photon counts, = 1 , and large photon counts, = l+ 10 ^(£’^/GeV)“^, which could potentially 
be achieved in a stacked halo analysis [22, 24]. The E~‘^ dependence is included as a crude model of 
the ganuna ray flux: lower energy gamma rays are more abundant. 

The preceding discussion yields a range of field strengths where measurements of cascade halo 
sizes are best suited to probing the IGMF. This range depends on the energy of the observed gamma 
rays, as illustrated in Fig. 11 , and taking the extremal values 1 GeV < E^ < 100 GeV leads to 
10“^^ G < Bq < 10“^^ G. However, the range also depends on photon fluxes. For instance, probing 
Bq > G requires a sufficient abundance of high energy photons with E^ > 10 GeV, which re¬ 

main in the MBC regime, and probing Bq < 10“^^ G requires a large number of low energy photons 
E^ < 10 GeV so as to beat down the effective angular resolution. More conservatively speaking, 
the range of magnetic field strengths that can be probed using cascade halo size measurements lies 
between 10“^^ — 10“^^ G and 10“^^ — 10“^^ G. This range is indicated on the parameter space plots 
of Fig. 12 as the blue shaded region. The gradient in the blue color indicates where the boundaries 
depend on the details of the analysis, just discussed. This range is independent of the magnetic field 
coherence length A as long as A is sufficiently large that the charged lepton probes an effectively 
homogeneous magnetic field, Le., X > De ^ 100 kpc. For smaller coherence length, the charged 
leptons do not follow helical arcs; instead their motion is that of a random walk, and our analysis 
breaks down. 

^This formula matches Fig. 2 of [43] for the case of on-orbit data (P6_V11) and front-converting events. 
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Figure 11. The reach of a gamma ray telescope like Fermi-LAT to probing an IGMF of strength Bq using 
cascade gamma rays of energy E^. We consider blazars located at both dg = I Gpc (blue) and 0.7 Gpc (red) 
from Earth. The two upper lines labeled 6 ^ De/R < I demarcate the boundary between the MBC regime, 
where a measurement of the halo size furnishes a measurement of the field strength, and the PH regime, where 
the halo size is insensitive to the field strength. The two middle curves labeled ©ext > ^Oqs indicate the 
field strength below which the halo size is smaller than the Fermi-LAT 68% confinement radius. Finally, we 
suppose that many cascade photons are observed = 1 10^{E^/GeV)~‘^, and the two lower curves, 

©ext > dOes / indicate the approximate field strength below which the halo cannot be resolved. 


Short Coherence Length Regime 

When the coherence length is small, information about spatial inhomogeneities is encoded in features 
of the halo at small angular scales. It is challenging to measure this halo substructure given limitations 
on angular resolution and photon flux. If one is only able to measure the large scale halo morphology, 
then gamma ray observations have a diminished capacity to probe magnetic helicity on small length 
scales. 

As a specific example, recall the study of helical magnetic fields in Case 4 of Sec. 4.4. In the left 
panel of Fig. 7, the halo of GeV gamma rays form a spiral around the blazar. For small A the angular 
scale of the spiral (separation between subsequent cycles) is smaller than the angular scale of the full 
halo. If we coarse grain on the scale of the halo, e.g. to model the finite detector resolution, the map 
would appear rotationally symmetric, but when the small scale behavior is resolved, the spiral can be 
seen. 

To study the small scale structure of the halo, we require not only good angular resolution, but 
also closely spaced gamma ray energies. Recall from Eqs. (2.3) and (2.5) that the mean free path of 
the TeV gamma ray is 

= (lo^) ' 

where E^ is the energy of the GeV gamma ray reaching Earth. Each TeV gamma ray emitted by the 
blazar can be viewed as a local probe of the magnetic field at the point where pair production occurs. 
If we want to probe the magnetic field on the length scale A, then we should employ gamma rays 
with energies Ei = E^ and E 2 — E^ + AE^ such that (Ei) — d^^ (£’ 2 ) ~ A. The optimal energy 
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separation is therefore 


AE^ ^ 1 GeV 
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(5.7) 


assuming AE^ <C E^ for which A {dds/dE^)AE^. Considerations of gamma ray flux prohibit 
one from taking AE^ arbitrarily small to probe arbitrarily small A. 

Since the Q-statistics are constructed from three different gamma ray energies, as in Eq. (3.23), 
the associated AE^ and E^ for a given statistic lead to a corresponding optimal coherence length 
Aopt via Eq. (5.7). The effects of inhomogeneities on smaller length scales, A < Aopt d^Q ( EJ ^) 

d^Q {E^ + AE^), are washed out. We have seen this behavior explicitly for Case 4 in Eig. 8 and Case 
5 in Eig. 10: as the coherence length decreases below A ^ Aopt, the Q-statistic stops growing and 
turns over or flattens out. To ensure that the optimal coherence length of a particular Q-statistic is 
small enough to probe magnetic field inhomogeneities on the scale A, we impose 


A > 10 Mpc 
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While a non-zero measurement of Q is of interest even for smaller A, it becomes difficult to draw 
a connection between the value of Q and the parameters of the underlying magnetic field. This is 
particularly evident in Eig. 8 where Q oscillates with varying A, and even sign[(5] does not neces¬ 
sarily correspond to the sign of the magnetic helicity. Presumably detailed modeling of the cascade 
development would be required to infer properties of the IGME at such small A. 

In the parameter space plots of Eig. 12 we indicate the coherence length below which oscilla¬ 
tions set in by writing “Q oscillations.” This boundary is drawn at the point where each of the Q- 
statistics has its first zero-crossing, determined numerically, but it is well-approximated by Eq. (5.8). 
Eor the energies and energy spacings that we have shown, the Q-statistics are best suited to probe 
coherence lengths in the range 10 — 100 Mpc. In order to probe shorter coherence lengths, Eq. (5.8) 
indicates that we need to consider Q-statistics built from larger E^ and smaller AE^. However, this 
limit is technically challenging: higher energy gamma rays are less abundant and this problem is 
exacerbated by small energy bins. 

Long Coherence Length Regime 

We have already seen in the analyses of Secs. 4.4 and 4.5 that the Q-statistic becomes small as 
the coherence length is made large. In this limit, the magnetic field is effectively homogeneous on 
the scale probed by the cascade, and information about the magnetic helicity, which is encoded in 
gradients, becomes inaccessible. In terms of the halo map, the characteristic “S”-like curve becomes 
flattened out (see the left panel of Eig. 7 or 9), and the transverse angular extent of the halo becomes 
small. As a consequence angular resolution is a challenge for a measurement of small Q since one 
needs not only to distinguish the halo from a point, but also distinguish the flattened S from a line. 

Recall that the triple product Qabc = x • fie can be written as in Eq. (3.24) 


Qabc = sin 'dab sin ipabc 


(5.9) 


where •dab is the angle between ha and h^, and ipabc is the angle between he and the plane spanned 
by the other two vectors. Because of the detector’s finite angular resolution, the unit vectors can only 
be localized on the sphere with a precision of S6qs. This translates into a comparable uncertainty on 
•dab and (fabc, and we estimate the precision with which the Q-statistic can be measured as 

SQabc ~ V COs2 -dab sin^ (fabc + sin^ 'dab COs2 (fabc ^(*68 . (5.10) 
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Figure 12. This sensitivity plot gives a rough indication of the region of magnetic field parameter space that 
can be probed with a gamma ray telescope like the Fermi-LAT using four different Q-statistics, Eq. (3.23). In 
the blue shaded regions, the magnetic field strength can be inferred from a measurement of the halo size. For 
smaller field strength the halo cannot be distinguished from a point source, and for larger field strength (PH 
regime), one may still be able to measure the halo size, but its connection to the field strength is obscured, 
see Fig. 4. In the red shaded region, the magnetic helicity can be inferred from a measurement of the halo 
shape via the Q-statistic. (These curves are defined with respect to the helical field configuration in Case 4 of 
Sec. 4.4, and we expect similar boundaries for more general helical configurations.) For larger A and smaller 
B the parity-odd signal is too small to measure given the assumed angular resolution, and for smaller A one 
may still be able to measure the halo shape via Q, but its connection to the coherence length is obscured, see 
Fig. 8. 
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The factor on the right side depends on the choice of gamma ray energies used to form the triplet as 
well as the specific magnetic field configuration being considered. As a rough estimate, we assume 
the small angle approximation and take i^ab ^ ^abc ^ ©ext(10 GeV) using Eq. (3.8). Thus we infer 
that the Fermi-LAT angular resolution is satisfactory to measure a Q-statistic satisfying 

Q>5Q^ 0ext(lO GeV) 69^8 ~ 3 x , (5.11) 

which is analogous to Eq. (5.4) for the halo size. As per the discussion surrounding Eq. (5.4), we 
expect that large photon counts can weaken the bound. Since this is a rough estimate, we consider 
6Q ^ 10“^, 10“^, and 10“^ in the following analysis. 

In Fig. 12 we also show the boundaries corresponding to the inequality in Eq. (5.11) for each 
of the energy combinations considered previously. The heavy red shaded region corresponds to Q > 
6Q ^ 10~^ where gamma ray observations of cascade halos should be able to measure the strength 
and helicity of the magnetic field. Moving outside of the red shaded region, the helicity measurement 
becomes more challenging, because either Q is too small to measure given an angular resolution 
comparable to the Fermi-LAT (large A regime) or Q has a complicated dependence on the magnetic 
helicity (small A regime). Outside of the blue shaded region, even the field strength measurement 
becomes challenging, because either the halo is too small to distinguish from a point source given 
the angular resolution (small Bq regime) or the halo size is insensitive to the field strength (large Bq 
regime). With better statistics one can presumably beat down the angular resolution issues to probe 
weaker fields. 

6 Summary and Discussion 

One day gamma ray observations may provide measurements of both the size and the shape of cascade 
halos. Whereas the halo size is tied to the average magnetic field strength, the halo shape can depend 
on other properties of the magnetic field, particularly its helicity. In this work we have attempted to 
illuminate the connection between magnetic helicity and halo shape by studying the development of 
the electromagnetic cascade for simplified magnetic field configurations. We have taken an analytic 
approach in which the three equations, (3.17), (3.18), and (3.19), are solved to find the trajectory of 
gamma rays reaching Earth. 

By acting as local probes of the magnetic field, blazar-induced cascade halos offer an interesting 
opportunity to measure properties of the IGMF. The physical parameters of the cascade and the 
detector resolution come together to determine the region of magnetic field parameter space that can 
be probed with this technique. By studying a particular helical magnetic field configuration with its 
wavevector oriented along the line of sight to the blazar, we have estimated the boundaries of this 
region in the two dimensional parameter space consisting of field strength Bq and coherence length 
A; see Fig. 12. The boundaries are partly determined by the anticipated detector resolution, and an 
improvement in angular resolution would allow access to smaller Bq and larger A. The boundaries 
are also depend on the measurement strategy, namely the energy combinations used to form the 
parity-odd Q-statistics via Eq. (3.23). For the four statistics considered here, we find that blazar- 
induced cascade halos are well-suited to probing helical magnetic fields with coherence length A > 
10 Mpc, see Eq. (5.8). Reaching smaller length scales would require more closely spaced energy 
combinations, and then one runs into issues with angular resolution and photon flux. 

There are a number of avenues for future work. First, we have considered only the simplest mag¬ 
netic field configurations in order to draw the connection between halo shape and magnetic helicity. 
These configurations are not likely to be realistic models of the true IGMF. It would be straight¬ 
forward to generalize our analysis and consider field configurations built from multiple modes in 
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Lateral Extent: 9 cos 0 [deg] 

Figure 13. A “smeared” halo map for Case 4 with parameters Bq = 10“^^ G, A = 700 Mpc, and = 1 Gpc. 
The same parameters correspond to the the green curve in Fig. 7 (left panel), which is also shown on this figure 
as a black curve. The colors denote different gamma ray energies: b < E^/ GeV < 10 (magenta), 10 < 
E^/ GeV < 15 (blue), 15 < E^/ GeV < 20 (green), 20 < E^/ GeV < 30 (yellow), 30 < E^/ GeV < 50 
(orange), and 50 < GeV < 100 (red). 


random orientations. In this case, it may become too difficult to solve Eqs. (3.17)-(3.19), and in¬ 
stead one would want to simulate the cascade using a shooting algorithm, such as in the analysis of 
Ref. [15]. 

Second, the parity-odd test statistic Q is not an automatic proxy for the magnetic helicity. We 
have seen in Sec. 4.4 that even the sign of Q won’t necessarily equal the sign of the magnetic helicity 
if the coherence length is small. (This was also pointed out in [37, 39] in the case of statistically 
homogeneous and isotropic magnetic fields.) The statistic Q can still be useful in this regime if the 
three energies at which it is evaluated are finely spaced, though then the photon counts may be small 
and the resulting error bars will be large. Alternately, now that we know qualitatively the effects of 
magnetic field helicity in the gamma ray pattern, it may be possible to devise improved statistics to 
detect these patterns. 

Third, our analysis neglects the stochastic nature of the cascade’s development. This simplifi¬ 
cation allows us to calculate the halo map i^{E^) that gives a deterministic connection between the 
energy of a gamma ray and its arrival direction on the sky. A more realistic model would account 
for the stochasticity in the gamma ray propagation distances, spectra, and so on. As a crude model 
of this effect, we calculate a halo map for Case 4 by replacing and in Eqs. (3.17)-(3.19) with 
rid^Q and r 2 D^ where ri and r 2 are randomly selected from a uniform probability distribution over 
the range 0.5 < ri, r 2 < 1.5. We calculate 100 such halo maps in this way for the same parameter 
set and show the results in Eig. 13. Although the degree of smearing is significant, the halo map 
retains the qualitative features that we have seen previously, namely it angular extent and parity-odd 
spiraling shape. 

Eourth, our analysis neglects the jet structure of the blazar, and assumes that the emission is 
isotropic. It is straightforward to include the effect of the jet in our semi-analytic formalism. Let 

^jet (^off 5 0off) ~ sin(0off )p(^off 5 0off) COs(0off ? 0off) (6.1) 
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be the orientation of the jet where 0off and (/)ofr give the polar and azimuthal offset angles from the 
line of sight. The orientation of the initial TeV gamma ray, v^, is given by Eq. (3.20). Then the jet 
criterion, 


arccos 


Vi 



^ ^jet 5 


( 6 . 2 ) 


ensures that the TeV gamma ray is emitted from within the cone of the jet. As a consequence of 
the jet criterion, only a portion of the full halo map becomes illuminated. This behavior is shown in 
Fig. 14 where we have reevaluated the halo maps from Case 4 using Eq. (6.2), and there is a similar 
effect for the other cases. When 0off > ^iet the Earth is not contained within the cone of the jet, an 
the blazar itself may not be seen. In this regime, the cascade photons would appear to contribute to 
the diffuse gamma ray flux. 
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Figure 14. The effect of the finite jet angle, ^jet, and offset from the line of sight, 6>off, where 0off = 0. The 
thin colored curves show the halo map for an isotropically emitting blazar, ^jet = 180°, and the thick curves 
show the halo map for a blazar jet with 6>jet = 5° and different offsets 6>off = 0°, 5°, and 8° in the three rows, 
respectively. For 6>off >10° none of the halo is visible with > 5 GeV. 
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